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Abstract —In this paper, we investigate power efficient resource 
allocation algorithm design for multiuser wireless communication 
systems employing a full-duplex (FD) radio base station for 
serving multiple half-duplex (HD) downlink and uplink users 
simultaneously. We propose a multi-objective optimization 
framework for achieving two conflicting yet desirable system 
design objectives, i.e., total downlink transmit power minimization 
and total uplink transmit power minimization, while guaranteeing 
the quality-of-service of all users. To this end, the weighted 
Tchebycheff method is adopted to formulate a multi-objective 
optimization problem (MOOP). Although the considered MOOP 
is non-convex, we solve it optimally by semidefinite programming 
relaxation. Simulation results not only unveil the trade-off between 
the total downlink and the total uplink transmit power, but also 
confirm that the proposed FD system provides substantial power 
savings over traditional HD systems. 

1. Introduction 

The exponential growth in high data rate communication 
has triggered a tremendous demand for radio resources such 
as bandwidth and energy. A relevant technique for reducing 
the energy and bandwidth consumption for given quality-of- 
service (QoS) requirements is multiple-input multiple-output 
(MIMO) as it offers extra degrees of freedom for more efficient 
resource allocation. However, the MIMO gain may be difficult 
to achieve in practice due to the high computational complexity 
of MIMO receivers m. As an alternative, multiuser MIMO 
(MU-MIMO) has been proposed as an effective technique for 
achieving the MIMO performance gain. In particular, in MU- 
MIMO, a transmitter with multiple antennas (e.g. a base station 
(BS)) serves multiple single-antenna users which shifts the 
computational complexity from the receivers to the transmitter 
m. Nevertheless, the spectral resource is still underutilized 
even if MU-MIMO is employed due to the traditional half¬ 
duplex (HD) operation of the BS, where uplink and downlink 
communication are separated in either time or frequency which 
leads to a significant loss in spectral efficiency. 

As a result, full-duplex (FD) wireless communications has 
recently received significant attention from both academia and 
industry due to its potential to double the spectral efficiency of 
the existing wireless communication systems O-lQ. In contrast 
to conventional HD transmission, FD enables simultaneous 
downlink and uplink transmission at the same frequency. 
However, in practice, a major challenge in FD communication 
is self-interference (SI) caused by the signal leakage from the 
downlink transmission to the uplink signal reception. Besides, 
the resource allocation for FD systems is complicated by the co¬ 
channel interference (CCI) caused by the uplink transmission to 
the downlink transmission. Thus, different resource allocation 
designs for FD systems were proposed and studied to overcome 
these challenges. The authors of (H investigated the end-to-end 
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outage probability of MIMO FD single-user relaying systems. 
In O, a resource allocation algorithm was proposed for the 
maximization of the end-to-end system data rate of multi¬ 
carrier MIMO FD relaying systems. In (61, massive MIMO 
was applied in FD relaying systems to facilitate SI suppression 
and to improve spectral efficiency. Simultaneous downlink and 
uplink transmission in small cells was studied in 0, where 
a suboptimal downlink beamformer was designed to improve 
the system throughput. However, the power consumption in FD 
systems has not been thoroughly studied in the literature so far 
( 1 - 0 . 

On the other hand, a large amount of work has been 
devoted to power efficient/green communications because of 
the pressure stemming from environmental concerns and the 
rapidly increasing cost of energy (Sl- lfTTl . In 0, downlink 
beamforming was proposed to minimize the total transmit 
power while guaranteeing secure communication for a single 
desired user. In m, both uplink and downlink beamforming 
were studied for total transmit power minimization in MU- 
MIMO systems. In (TOl . a power efficient scheduling scheme 
was designed to minimize the power consumption of multihop 
wireless links. In im, the maximization of the energy efficiency 
of multi-carrier systems was studied. However, all of the above 
works focus on HD systems and the obtained results may 
not be applicable to FD communication systems. In particular, 
total downlink and total uplink transmit power minimization 
are conflicting design objectives in FD systems. However, 
the intricate trade-off between downlink and uplink power 
consumption in FD systems has not been studied in the literature 
ii-illll. Such a trade-off does not exist in HD systems since 
downlink and uplink transmission are separated either in the 
time domain or in the frequency domain. Besides, the single¬ 
objective optimization frameworks in (Sl- lfTll are not suitable 
for studying the trade-off between downlink and uplink power 
consumption in FD systems. 

Motivated by the aforementioned observations, we propose 
a multi-objective optimization problem (MOOP) formulation 
which minimizes the total downlink and the total uplink transmit 
powers subject to QoS constraints for MU-MIMO FD wireless 
communication systems. The proposed MOOP is formulated 
by adopting the weighted Tchebycheff method and solved 
optimally by semidefinite programming (SDP) relaxation. 

H. System Model 

In this section, we present the adopted MU-MIMO channel 
model for FD wireless communication systems. 

A. Notation 

We use boldface capital and lower case letters to denote 
matrices and vectors, respectively. A^, Tr(A), and Rank(A) 
represent the Hermitian transpose, trace, and rank of matrix A, 


respectively; and represent the inverse and Moore- 

Penrose pseudoinverse of matrix A, respectively; A ^ 0 
indicates that A is a positive semidefinite matrix; I at is the 
N X N identity matrix; denotes the set of all A x M 

matrices with complex entries; denotes the set of all A x A 
Hermitian matrices; 1*1 and H-H denote the absolute value of a 
complex scalar and the Euclidean vector norm, respectively; 
£{■} denotes statistical expectation; the circularly symmetric 
complex Gaussian distribution with mean fi and variance is 
denoted by and ^ stands for “distributed as”. 



B. Multiuser System Model 

We consider a multiuser communication system. The system 
consists of an FD radio BS, K downlink users, and J uplink 
users, cf. Figure [T] The FD radio BS is equipped with At > 
1 antennas for simultaneous downlink transmission and uplink 
reception in the same frequency bancfl. The A + J users are 
single-antenna HD mobile communication devices to ensure low 
hardware complexity. The downlink and the uplink users are 
scheduled for simultaneous uplink and downlink transmission. 
Besides, we assume that the global channel state information 
(CSI) of all users is perfectly known at the BS for resource 
allocation. 

C. Channel Model 

We focus on a frequency flat fading channel. The received 
signals at downlink user k G {1,..., A} and the FD radio BS 
are given by 

K 

Vk^ = hfxfe + hf xm 
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self-interference 

respectively, where x/. G denotes the transmit signal 

vector from the FD radio BS to downlink user k. The downlink 
channel between the FD radio BS and user k is denoted by G 
C^txi ^ ^ ^ represents the channel between uplink user 

j and downlink user k. Variables d^^ and Pj are the transmit 
data and the transmit power sent from uplink user j to the FD 
radio BS, respectively. Vector gj G denotes the channel 

between uplink user j and the FD radio BS. Matrix Hsi G 
£NtxNt denotes the SI channel at the FD radio BS. Variables 
h/c, fj,k^ gjs and Hsi capture the joint effect of path loss and 
small scale fading, z ^ CA(0, and ^ CA(0, 

represent the additive white Gaussian noise (AWGN) at the FD 
radio BS and downlink user k, respectively. In ([T]), the term 
represents the multiuser interference (MUI) in 
the downlink channel, and the term \/^fj,kd^^ denotes 


Fig. 1. A multiuser communication system with a full-duplex (FD) radio 
base station (BS), K = 1 half-duplex (HD) downlink users, and J = 1 HD 
uplink users. The BS is equipped with A^t antennas for simultaneous uplink 
and downlink communication. 


the CCI created ^ the uplink users to downlink user k. In (O, 
the term Hsi X]/c=i represents the SI. 

In each scheduling time slot, the FD radio BS transmits 
A independent signal streams simultaneously at the same 
frequency to the A downlink users. In particular, the transmit 
signal vector 

Xfe = (3) 

is adopted at the BS, where d^^ G C and Wk G are 

the information bearing signal for downlink user k and the 
corresponding beamforming vector, respectively. Without loss 
of generality, we assume £{\d^^\‘^} = £{\d^^\‘^} = l,V/c G 
{l,...,A},jG{l,...,J}. 

III. Problem Formulation 

In this section, we first define the adopted performance metric 
for the considered multiuser communication system. Then, we 
formulate the downlink and uplink power allocation problems. 
The number of antennas at the FD radio BS is assumed to be 
larger than the number of uplink users to facilitate uplink signal 
detection, i.e.. At > J. 


A. Performance Metrics 

The receive signal-to-interference-plus-noise ratio (SINR) at 
downlink user k is given by 


pDL 

J- k 
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E |hfw„^ 
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On the other hand, the receive SINR of uplink user j at the 
FD radio BS is given by 


pUL ^ 


A Is 






K 


(5) 


■ X] |vyHsiWfc|2 

k=l 


where e is the receive beamforming vector for 

decoding of the information received from uplink user j. In this 
paper, zero-forcing receive beamforming (ZF-BF) is adopted. 
We note that ZF-BF closely approaches the performance of 
optimal minimum mean square error beamforming (MMSE-BF) 
when the noise term is not dominating eB or the number 
of antennas is sufficiently large ||6l. Besides, ZF-BF facilitates 
the design of a computational efficient resource allocation 


^We note that circulator based FD radio prototypes which can transmit and 
receive signals simultaneously on the same antennas have been demonstrated 
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^We note that the noise power at the BS is not expected to be the dominating 
factor for the system performance since the BS is usually equipped with a high 
quality low-noise amplifier (LNA). 

















algorithm. Hence, the receive beamformer for uplink user j is 
chosen as 


Vi = (uiQ^)", 


( 6 ) 


where uj 



(i-i) (J-j) 

Q = [gi,---,gj]- 


(Q^Q)-iQ^, and 


B. Optimization Problem Formulation 

We first study the problem formulation of two desirable 
system design objectives for the considered FD wireless system. 
Then, we investigate the two system objectives jointly under the 
framework of multi-objective optimization. The first considered 
problem is designed to minimize the total downlink transmit 
power and given by 

Problem 1 (Total Lk?wnlink Transmit Power Minimization): 


minimize 

^k,Pj 


El 

k=l 


W/c 


s.t. Cl: 


C2: 




K 

E 

m^k 


hfw^ 


i=i 


> 

- reqfc > 


Vfc, 




-Pilgfvi 


5^^ + 0-2 


UL 
2 - - req^. 


> r 


Vj, 


C3: Pj > 0, Vj, 
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where we define = ^Pr\gr^j\‘^ + Z] for 

r^j k=l 

notational simplicity. Besides, > 0 and rUL > 0 are the 
minimum required SINRs for downlink user k G { 1 ,..., Ff} 
and uplink user j G respectively. 

Constraint C3 in © is the non-negative power constraint 
for uplink user j. Problem 1 is a non-convex problem due to 
constraints Cl and C2. We note that Problem 1 tries to minimize 
only the total downlink transmit power under constraints C1-C3 
without regard for the consumed uplink transmit powers. 

The second desirable system objective is the minimization 
of total uplink transmit power and can be mathematically 
formulated as 

Problem 2 (Total Uplink Transnjit Power Minimization): 


minimize 

^k,Pj 




s.t. Cl - C3. 


( 8 ) 


resort to multi-objective optimization. In the literature, multi¬ 
objective optimization is often adopted to study the trade-off 
between conflicting system design objectives via the concept 
of Pareto optimality. A point is Pareto optimal if there is no 
other point that improves at least one of the objectives without 
detriment to another objectively. In order to capture the complete 
Pareto optimal set, we formulate a third optimization problem 
to investigate the trade-off between Problem 1 and Problem 2 
by using the weighted Tchebycheff method □a. The resulting 
problem formulation is given as: 

Problem 3 (Multi-Objective Optimization): 

minimize max < (Qi{^k,Pj)-Q*)] 


s.t. Cl - C3, 


(9) 


where Qi{wk,Pj) = Ef=il|wfe||^ and Q 2 {yfk,Pj) = 
Constant Q* is the optimal objective value of the 
i-th problem. In particular, variable Q* is assumed to be a 
constant for resource allocation. Variable > 0, E.A, = 1 , 
specifies the priority of the i-th objective compared to the other 
objectives and reflects the preference of the system operator. By 
varying A^, we can obtain a complete Pareto optimal set which 
corresponds to a set of resource allocation policies. Thus, the 
operator can select a proper resource allocation policy from 
the set of available policies. Compared to other formulation 
methods for handing MOOPs in the literature (e.g. the weighted 
product method and the exponentially weighted criterion mi 
the weighted Tchebycheff method can achieve the complete 
Pareto optimal set with a lower computational complexity, 
despite the non-convexity (if any) of the considered problem. 
It is noticed that Problem 3 is equivalent to Problem i when 
A^ = 1 and Xj = 0, Vi 7 ^ j. Here, equivalent means that both 
problem formulations have the same optimal solution. 


IV. Solution of the Optimization Problems 

The optimization problems in (|7]), ([ 8 ]), and (|9]) are non-convex 
problems due to the non-convexity of constraints Cl and C2. In 
order to solve the problems efficiently, we recast Problems 1, 
2, and 3 as convex optimization problems via SDP relaxation. 

A. Semidefinite Programming Relaxation 

For facilitating the SDP relaxation, we define 

Wfc = Wfcwf , Hfc = hfehf , Gj = gjgf, Vj = Vjvf, (10) 


Problem 2 targets only the minimization of the total uplink 
transmit power under constraints C1-C3 without regard for the 
consumed downlink transmit power. 

The two above system design objectives are desirable for the 
system operator and the users, respectively. However, in FD 
wireless communication systems, these objectives conflict with 
each other. On the one hand, the uplink transmission causes CCI 
to the downlink transmission. Hence, the FD radio BS should 
transmit with a high power to ensure the QoS requirements of 
the downlink users. On the other hand, a high downlink transmit 
power results in strong SI to the uplink reception. Hence, the 
uplink users should transmit with a high power to satisfy the 
minimum required receive SINRs of the uplink users at the FD 
radio BS. However, this would in turn cause high CCI and give 
rise to successive increases in transmit power for both uplink 
and downlink transmission. To overcome this problem, we 


and rewrite Problems 1-3 into the following equivalent 
formulations: 

Transformed Problem 1: ^ 

minimize Tr (W/c) 


s,. + vt, 

C2: > ly^ + 0-2 Tr(Vj), Vj, 


pUL 
^ req. 
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C3: Pj > 0, Vj, 

C4: W/c E 0, \fk, 

C5: Rank(W/c) < 1, \fk, (11) 


^Please refer to □a for a detailed discussion of Pareto optimality. 









where Wfe ^ 0, Wk € and Rank(Wfe) < 1 in 

dnii are imposed to guarantee that W/c = w/cw|^ holds after 
optimization, and 

K J 

= Y. Tv(HkWm) + and ( 12 ) 

m^k i=l 

J K 

= E^rTr(V,G,)+^Tr(WfcH"V,Hsi). (13) 

r^j k=l 

Transformed Problem 2: j 

minimize > Pj 

Wfce]H^T,p^. ^ 

s.t. Cl - C5. (14) 

Transformed Problem 3: 

minimize t 

WkEM^T^Pj^t 

S.t. Cl - C5, 

C 6 : A,(Q,-Q*)<t,ViG{l,2}, (15) 

where t is an auxiliary optimization variable and ilSi is the 
epigraph representation of (O. 

Since Problem 3 is a generalization of Problems 1 and 2, 
we focus on the solution of Problem 3. Transformed Probl^ 
3 is a non-convex problem due to the rank-one constraint C5 
and solving such a rank-constrained problem is known to be 
NP-hard d. Hence, we relax constraint C5: Rank(W/c) < 1 
by removing it from the problem formulation, such that the 
considered problem becomes a convex SDP given by 

minimize t 

s.t. C1-C4,C6. (16) 

The relaxed convex problem in ([16]) can be solved efficiently by 
standard convex program solvers such as CVX ca. Besides, 
if the solution obtained for a relaxed SDP problem is a rank- 
one matrix, i.e., Rank(W/c) = 1 for W/^ 7 ^ 0 , \fk, then it is 
also the optimal solution of the original problem. Otherwise, in 
general, the solution of ([T^ serves as a lower bound for the 
original problem (fTSl) due to the larger feasible solution set in 

(HSI). 

Next, we reveal a sufficient condition for obtaining a rank-one 
solution W k for the relaxed version of Transformed Problem 3 
in the following theorem. 

Theorem 1: Under the condition that the channel vectors 
of the downlink users, h/c,/c G and the uplink 

users, gj,j G J}, can be modeled as statistically 

independent random variables, the solution of ([T^ is rank-one, 
i.e., Rank(W/c) = 1 for W/c 7 ^ 0 , V/c, with probability one. 

Proof: Please refer to the Appendix. ■ 

In other words, whenever the channels satisfy the condition 
stated in Theorem 1, the optimal beamformer of ([9]) can be 
obtained with probability one. 

Remark 1: Exploiting Theorem 1, Problem 1 and Problem 2 
can also be solved optimally via SDP relaxation. In particular, 
for Ai = 1 , A 2 = 0 , and Q\ = 0 , the problem in ([16]) becomes 
Problem 1; for A 2 = l,Ai = 0, and Q 2 = 0, the problem 
in ([16]) becomes Problem 2. We note that setting (5* = 0 in 
Problem 3 is used only for recovering the solution of Problems 
1 and 2. However, this does not imply that the optimal value 
of Problem i, Vi G {1, 2}, is equal to 0. 


TABLE I 

System parameters. 


Carrier center frequency 

1.9 GHz 

System bandwidth 

200 KHz 

Path loss exponent 

3.6 

Reference distance 

30 m 

SI cancellation 

-80 dB |3l 

Downlink user noise power, 

—83 dBm 

Uplink BS noise power, 

-110 dBml 

BS antenna gain 

10 dBi 


V. Results 

In this section, we investigate the performance of the 
proposed multi-objective optimization based resource allocation 
scheme through simulations. The simulation parameters are 
specified in Table jl] There are K = 3 downlink users and J = 8 
uplink users in the cell. The users are randomly and uniformly 
distributed between the reference distance and the maximum 
service distance of 250 meters. The FD radio BS is located at 
the center of the cell and equipped with antennas. The small 
scale fading of the downlink channels, uplink channels, and 
inter-user channels is modeled as independent and identically 
distributed Rayleigh fading. The multipath fading coefficients 
of the SI channel are generated as independent and identically 
distributed Rician random variables with Rician factor 5 dB. 

A. Transmit Power Trade-off Region 

In Figure [2] we study the trade-off between the downlink 
and the uplink total transmit powers for different numbers of 
antennas at the FD radio BS. The trade-off region is obtained 
by solving ([16]) for different 0 < A^ < l,Vi G {1,2}, i.e., the 
A^ are varied uniformly using a step size of 0.01 subject to 
Ai = 1. Without loss of generality, we assume all downlink 
users and all uplink users require the same SINRs, respectively, 
i.e., = rDL = 10 dB and = 6 dB. It can 

be observed from Figure 2 that the total uplink transmit power 
is a monotonically decreasing function with respect to the total 
downlink transmit power. In other words, minimizing the total 
uplink power consumption leads to a higher power consumption 
in the downlink and vice versa. This result confirms that the 
minimization of the total uplink transmit power and the total 
downlink transmit power are conflicting objectives. For the case 
of Nt = 10, 10.9 dB in uplink transmit power can be saved 
by increasing the total downlink transmit power by 6.5 dB. 
In addition. Figure [2] also indicates that a significant amount 
of power can be saved in the FD system by increasing the 
number of antennas. This is due to the fact that the extra 
degrees freedom offered by the additional antennas facilitate 
a more power efficient resource allocation. However, there is 
a diminishing return in the power saving as the numbers of 
antennas at the FD BS increase due to channel hardening ina. 

For comparison, we also consider a baseline multiuser system 
with a HD radio BS. In particular, the HD radio BS is 
equipped with Np antennas which are utilized for transmitting 
or receiving in non-overlapping equal-length time intervals. As 
a result, both SI and CCI are avoided. In order to have a 
fair comparison between the FD and the HD systems, we set 

log2(l+r?4) = l/21og2(l + rP 4 -™) andlog2(l+ry4) = 

^The downlink user noise power is higher than the uplink BS noise power 
since, unlike cheap mobile devices, the BS is usually equipped with a high 
quality LNA and a high resolution analog-to-digital converter (ADC) E). 
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Fig. 2. Average system objective trade-off region achieved by the proposed 
optimal resource allocation scheme. The double-sided arrows indicate the power 
saving due to additional antennas. 


Fig. 3. Average power consumption (dBm) versus the minimum required 
downlink SINRs (dB) for FD and HD radio BS. 


1/2 log 2 (l + rJ^q. ™). Thus, for the HD radio BS, the required 
SINRs for the downlink and uplink users are = 

(1 + rfekrUL-HD = (1 + pUL )2 _ respectively. 
In addition, the power consumption of downlink and uplink 
transmission for the HD system is divided by two since either 
uplink or downlink transmission is performed at a given time. 
Then, we optimize W/c to minimize the downlink total transmit 
power, and we optimize Pj to minimize the uplink total transmit 
power for the optimal MMSE receiver at the HD radio BS 
na. The power consumption for the baseline HD system is 
represented by a single point in Figure [2| for a given number of 
antennas. It can be observed that the power consumption of the 
baseline system lies above the curve of the corresponding FD 
system. In the regions which are spanned by the dashed lines 
and the corresponding FD system performance curves, the FD 
systems are more power efficient than the corresponding HD 
systems in both uplink and downlink. In particular, for the case 
of Nt = 10, the FD system facilitates more than 6 dB power 
saving for downlink transmission and more than 5 dB power 
saving for uplink transmission compared to the corresponding 
HD system. 


B. Average Transmit Power 

In Figure [3l we investigate the power consumption of 
downlink and uplink users for identical minimum required 
downlink SINRs, T^q = Tj^q, and identical minimum required 
uplink SINRs of = 6 dB. The FD radio BS and the 

baseline HD radio BS are both equipped with = 10 
antennas. We select the resource allocation policy with Ai = 0.1 
and A 2 = 0.9 which indicates that the system operator attaches 
a higher priority to total uplink transmit power minimization. 
It can be observed that both the downlink and the uplink 
power consumption increase with Tj^q. However, the downlink 
power consumption grows more rapidly than the uplink power 
consumption. The reason behind this is twofold. First, a higher 
downlink transmit power is required to fulfill more stringent 
downlink QoS requirements. Second, the SI becomes more 
severe for higher downlink transmit powers. Specifically, since 
the FD radio BS tries to keep the total uplink transmit power 
low to avoid strong CCI for increasing downlink transmit power. 


more degrees of freedom at the FD radio BS have to be utilized 
to suppress the SI to improve the receive SINR of the uplink 
users. As a result, fewer degrees of freedom are available for 
reducing the downlink transmit power consumption causing 
even higher downlink transmit powers. Besides, Figure 3 also 
confirms that both downlink and uplink transmission in the 
proposed FD system are more power efficient compared to the 
HD system. 


VI. Conclusions 

In this paper, we studied the power allocation algorithm 
design for MU-MIMO wireless communication systems with 
an FD BS. The algorithm design was formulated as a non- 
convex MOOP employing the weighted Tchebycheff method. 
The proposed MOOP was solved optimally by SDP relaxation. 
Simulation results not only unveiled the trade-off between the 
total downlink and the total uplink transmit powers, but also 
confirmed that the proposed FD system enables substantial 
power savings compared to a HD system. 


Appendix - Proof of Theorem 1 

The relaxed version of Transformed Problem 3 in (O is 
jointly convex with respect to the optimization variables and 
satisfies the Slater’s constraint qualification. Therefore, strong 
duality holds and solving the dual problem is equivalent to 
solving the primal problem ini- For obtaining the dual problem, 
we first need the Fagrangian function of the primal problem in 
(O which is given by 

c(Wk,Pj,t,^k,'>Pj,Yk,aj,0i^ (17) 

= ^IV(RfeW,)-^Iv(Wfe(Y, + |^))+A, 

k=l k=l 
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^ = t{\-Y,Qi)-e^\iQi + e2\2{Y,Pj-Ql)^ (19) 

i=l i=l 

K J 

Kk = ^ ^rnH^ + + 6>iAilArT- (20) 

m^k j=l 

Here, A denotes the collection of terms that only involve 
variables that are independent of W/^. ^k, OLj, and Oi are 
^ Lagr^ge multipliers associated with constraints Cl, C2, 
C3, and C 6 . Matrix Y/. G is the L^ange multiplier 

matrix for the positive semidefinite constraint C4 on matrix W/c. 
Thus, the dual problem for the SDP relaxed problem in (O is 
given by 

maximize minimize £(W/e, P.-, t, £/e, Y/., a.-, 

^fc,'0z,a>O WfcGlHI^T V / 

0i,^k^O p-,t 

2 

s.t. Y.^i = l. (21) 

i=l 

The constraint = 1 is imposed to guarantee a bounded 

solution of the dual problem ifTTll . 

We reveal the structure of the optimal W/^ of (O by studying 
the Karush-Kuhn-Tucker (KKT) conditions. For convenience, 
we define ^ {W*,P;,r} and ^ Y*,a*,^*} 

as the sets of the optimal primal and dual variables of ( 121 ]) . 
In the following, we focus on the KKT conditions which are 
needed for the proof: 

Y^PO, ( 22 ) 

Y*W^ = 0 , and (23) 

(24) 

J- req^ 

where is obtained by substituting the optimal dual variables 
into (l20l) . By exploiting (|23]), we have Y^W^ = 0 ^ 
Rank(Y^) + Rank(W^) < At. By exploiting the fact that 
Rank(W^) > 1 for 7 ^ 0 , we have 

Rank(Y*) < At - 1. (25) 

If Rank(Y^) = Nt — 1, then the optimal beamforming matrix 
7 ^ 0 must be a rank-one matrix. In order to further reveal 
the structure of Y^, we first show by contradiction that is a 
positive-definite matrix with probability one under the condition 
stated in Theorem 1. Let us focus on the dual problem in (121]) . 
For a given set of optimal dual variables, $*, the optimal uplink 
transmit power, P*, and the optimal auxiliary variable t*, the 
dual problem in ([ 2 T]) can be written as 

minimize £ (w,, P/, P, , Y*, a*, . ( 26 ) 

Suppose R^ is not positive-definite. In this case, we can choose 
W/c = rw/cW^ as one of the optimal solution of ([26b . where 
r > 0 is a scaling parameter and w/c is the eigenvector 
corresponding to a non-positive eigenvalue p/c < 0 of R^, 
i.e., R^w/c = p/cW/c. We substitute W/c = rwkW^ and 
Rlwk = Pk^k into (l26b which leads to 

^ lV(pferw,wf) - r ^ IV (w^wf (Y^ + |^)) + A. 

k=l k=l 

(27) 

In (l27b . the first term, X]/c=i Tr(p/crw/cw|^) is not positive. 
As for the second term, since the channel vectors of gj and 


h/c are assumed to be statistically independent, and from (l 22 b . 
we obtain ^ Tr fv^/c^vf (Y^ -f 1^^)) > 0. By setting 

r ^ 00 , the term -r Tr ^v^/c^vf (Y^ + j ^ - 00 . 

In this case, the dual optimal value becomes unbounded from 
below. However, the optimal value of the primal problem ([T 6 b is 
non-negative due to constraint C 6 . Thus, strong duality cannot 
hold which leads to a contradiction. Therefore, R^ is a positive- 
definite matrix with probability one, i.e., Rank(R^) = At. By 
applying (l24b and a basic inequality for the rank of matrices, 
we have 

Rank(Y^) + Rank(C^^) (28) 

J- req^ 

> Rank {Yl + Ck^) = Rank(RJ) = Yt 

J- req^ 

^ Rank(Y^) > At - 1. 

Combining the inequalities in (l25b and ([28b, we obtain that 
Rank(Y^) = — 1. Thus, Rank(W^) = 1 holds with 

probability one. 

References 

[1] J. Zhang, M. Matthaiou, G. K. Karagiannidis, H. Wang, and Z. Tan, 
“Gallager’s Exponent Analysis of STBC MIMO Systems over 77-/1 and 
K,-fi Fading Channels,” IEEE Trans. Commun., voL 61, no. 3, pp. 1028- 
1039, Mar. 2013. 

[2] D. Gesbert, M. Kountouris, R. Heath, C.-B. Chae, and T. Salzer, “Shifting 
the MIMO Paradigm,” IEEE Signal Process. Mag., vol. 24, no. 5, pp. 36- 
46, Sep. 2007. 

[3] D. Bharadia, E. McMilin, and S. Katti, “Full Duplex Radios,” in Proc. 
ACM SIGCOMM, Aug. 2013, pp. 375-386. 

[4] H. Suraweera, I. Krikidis, G. Zheng, C. Yuen, and P. Smith, “Low- 
Complexity End-to-End Performance Optimization in MIMO Full-Duplex 
Relay Systems,” IEEE Trans. Wireless Commun., vol. 13, no. 2, pp. 913- 
927, Feb. 2014. 

[5] D. W. K. Ng, E. S. Lo, and R. Schober, “Dynamic Resource Allocation in 
MIMO-OFDMA Systems with Full-Duplex and Hybrid Relaying,” IEEE 
Trans. Commun., vol. 60, no. 5, pp. 1291-1304, May 2012. 

[ 6 ] H. Q. Ngo, H. Suraweera, M. Matthaiou, and E. Larsson, “Multipair Full- 
Duplex Relaying with Massive Arrays and Linear Processing,” IEEE J. 
Select. Areas Commun., vol. 32, no. 9, pp. 1721-1737, Sep. 2014. 

[7] D. Nguyen, L.-N. Tran, P. Pirinen, and M. Latva-aho, “On the Spectral 
Efficiency of Full-Duplex Small Cell Wireless Systems,” IEEE Trans. 
Wireless Commun., vol. 13, no. 9, pp. 4896-4910, Sep. 2014. 

[ 8 ] D. W. K. Ng, R. Schober, and H. Alnuweiri, “Power Efficient MISO 
Beamforming for Secure Layered Transmission,” in Proc. IEEE Wireless 
Commun. and Networking Conf., Apr. 2014, pp. 422-427. 

[9] M. Schubert and H. Boche, “Iterative Multiuser Uplink and Downlink 
Beamforming under SINR Constraints,” IEEE Trans. Signal Process., 
vol. 53, no. 7, pp. 2324-2334, Jul. 2005. 

[10] A. Lau and S. Cui, “Joint Power Minimization in Wireless Relay 
Channels,” IEEE Trans. Wireless Commun., vol. 6 , no. 8 , pp. 2820-2824, 
Aug. 2007. 

[11] D. W. K. Ng, E. S. Lo, and R. Schober, “Energy-Efficient Resource 
Allocation in OFDM A Systems with Large Numbers of Base Station 
Antennas,” IEEE Trans. Wireless Commun., vol. 11, no. 9, pp. 3292-3304, 
Jul. 2012. 

[12] D. Tse and P. Viswanath, Eundamentals of Wireless Communication, 
1st ed. Cambridge University Pres, 2005. 

[13] R. T. Marler and J. S. Arora, “Survey of Multi-objective Optimization 
Methods for Engineering,” Structural and Multidisciplinary Optimization, 
vol. 26, pp. 369-395, Apr. 2004. 

[14] B. GMner and J. Matousek, Approximation Algorithms and Semidefinite 
Programming. Springer Science & Business Media, 2012. 

[15] M. Grant and S. Boyd, “CVX: Matlab Software for Disciplined Convex 
Programming, version 2.1,” [Online] http://cvxr.com/cvx Mar. 2014. 

[16] H. Holma and A. Toskala, LTEfor UMTS-OFDMA and SC-EDMA Based 
Radio Access. John Wiley & Sons, 2009. 

[17] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge 
University Press, 2004. 








